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Abstract 

We consider a scalar-tensor model of dark energy with kinetic and Gauss Bon- 
net couplings. We study the conditions for the existence of quintessential and 
phantom power-law expansion, and also analyze these conditions in absence of 
potential (closely related to string theory). A mechanism to avoid the Big Rip 
singularity in various asymptotic limits of the model has been studied. It was 
found that the kinetic and Gauss-Bonnet couplings might prevent the Big Rip 
singularity in a phantom scenario. The autonomous system for the model has 
been used to study the stability properties of the power-law solution, and the 
centre manifold analysis was used to treat zero eigenvalues. 

PACS 98.80.-k, 95.36+x, 04.50.kd 



1 Introduction 

According to recent observations the current universe undergoes a phase of accel- 
erated expansion, due to domination of dark energy (DE) over the matter content 
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of the universe (usual barionic and dark matter) [TJ [2J, |3j HJ El E] . The astrophysical 
observations indicate that the equation of state parameter for dark energy lies in a 
narrow region around w — — 1, which might include values smaller than —1. The 
current observational data are in agreement with the simplest possibility of the cos- 
mological constant as the source of DE, but there is no mechanism to explain its 
smallness (expressed in Planck units) in contradiction with the expected value as the 
vacuum energy in particle physics [7], [E]; and observational data also show a better 
fit for a redshift dependent equation of state. Despite the variety of DE models, it is 
however difficult to fulfill all observational requirements like the observed value of the 
equation of state parameter (EoS) of DE, w ~ —1, the current content of DE relative 
to that of dark matter (known as coincidence problem), and the estimated redshift 
transition between decelerated and accelerated phases, among others. Among the 
models used to explain the DE (for review see [3 El O [101 El Q2]), the scalar-tensor 
theories are some of the most studied, not only because they contain direct couplings 
of the scalar field to the curvature, many of them predicted by fundamental theories 
like Kaluza-Klein and string theories [13], [TJ], but also because the scalar-tensor 
models fulfill in principle many of the above requirements. 

In the present work we consider a string and higher-dimensional gravity inspired 
scalar-tensor model, with non minimal kinetic and Gauss Bonnet (GB) couplings, to 
study late time cosmological dynamics and board some issues of dark energy. These 
terms are present in the next to leading a' corrections in the string effective action 
(where the coupling coefficients are functions of the scalar field) [T5] , [T6] and have 
the notorious advantage that lead to second order differential equations, preserving 
the theory ghost free. 

Some late time cosmological aspects of scalar field model with derivative couplings 
to curvature have been considered in [17], [321 [332 ED] , [21], [22]- On the other hand, 
the GB invariant coupled to scalar field has been extensively studied. In [23] the GB 
correction was proposed to study the dynamics of dark energy, where it was found 
that quintessence or phantom phase may occur in the late time universe. Different 
aspects of accelerating cosmologies with GB correction have been also discussed in 
[2~i] . [25] . [26] . [27], [28], [29]. The modified GB theory applied to dark energy has 
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been suggested in [30], and further studies of the modified GB model applied to late 
time acceleration, have been considered among others, in [31], [32], [33], [31], [33] . 
The combined effect of GB and kinetic coupling to curvature in the context of dark 
energy, has been considered in [3SJ ETJ EH ES] . In [3S] solutions with Big Rip and 
Little Rip singularities have been considered, in [37] 138] the reconstruction of different 
cosmological scenarios, including known phenomenological models has been studied, 
and in [39] some exact solutions have been found. 

In the present paper we focus on the power-law solutions of the quintessence and 
phantom types, in the case of late-time cosmology with scalar field dominance. We 
analyze the conditions to avoid the Big Rip singularity presented in phantom power- 
law solutions. The autonomous system for the model has been considered to find the 
restrictions on the parameter space of the model satisfying the conditions of stability. 
In section II we introduce the model and give the general equations, which are then 
expanded on the FRW metric to study the different power-law solutions. In section 
III we analyze the conditions to evade the Big Rip singularity in different scenarios. 
In section IV we introduce the dynamical variables and analyze the stability of the 
power-law solutions. Concluding remarks are given in section V. 



2 The model and power-law solutions 

We consider the following action that contains the Gauss Bonnet coupling to the 
scalar field and kinetic couplings to curvature (such terms are present in the leading 
a' correction to the string effective action [16J). 



S = d 4 xy/^g 



IQuG 2 

(2.1) 

where 7 = ±1 (+1 for the standard scalar field and —1 for the phantom scalar), 
G^y = R^ v — ^g^R, Q is the 4-dimensional GB invariant Q = R 2 — AR^R^ + 
R l ivp<jR ilvp0 ' ■ The coupling Fi(4>) has dimension of (length) 2 , and the coupling F 2 (0) 
is dimensionless. Note that compared with the more general action that leads to the 
second-order equations of motion (in metric and scalar field) [ID] , we are neglecting 
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derivative terms that are not directly coupled to curvature, of the form □0<9 A1 0<9 M 
and (d ll <j)d fl (j)) 2 , which is is acceptable in a cosmological scenario with accelerated 
expansion. The properties of the GB invariant guarantee the absence of ghost terms 
in the theory. Hence, the equations derived from this action contain only second 
derivatives of the metric and the scalar field. 



By varying Eq. (2.1 ) with respect to metric we derive the gravitational field equations 
given by the expressions 

i<W - = {Tr + Tfi>) (2.2) 

where k 2 = 8nG, T™ is the usual energy-momentum tensor for the matter component, 
the tensor T^ u represents the variation of the terms which depend on the scalar field 
4> and can be written as 

Tfiu = T*, + T% + T% B (2.3) 

where T^ v corresponds to the variations of the standard minimally coupled terms, 
T^ v comes from the kinetic coupling, and T^ U B comes from the variation of the cou- 
pling with GB. Due to the kinetic coupling with curvature and the GB coupling, the 
quantities derived from this energy-momentum tensors will be considered as effective 



ones. The respective components of the energy-momentum tensor (2.3) are given by 
It = 'WV^ - i 7 ^V A 0V A - 9flu V(<P) (2.4) 

T% = U pv - \g, v R\ F!(0)V A 0V A + ^,V A V A (F x (0)V 7 0V^) 



- ^(V M V, + V„V M ) (F!(0)V A 0V A 0) + RF x {<t>)V„cjN v cj) 

- 2F 1 (cf>) + J R, A V A 0V At 0) + ^i? A7 F 1 (0)V A 0V 
+ V A V^ (F!(0)V A 0V,0) + V A V, (F!(0)V A 0V M 0) 

- V A V A (F^V^V^) - <?^V A V 7 (F!(0)V A 0V^) 



(2.5) 



and 



T° B =4(iV,V u F 2 (<p)}R - gia ,\y p V p F 2 {<f>)]R - 2[VV/ 2 (^)]^ - 2[V p V u F 2 (<f))}R up 
+ 2[V p WF 2 {<f>)]R ltv + 2g, u [V^F 2 (<p)]R pa - 2[VV ff F 2 (^)] v) 

(2.6) 
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In this last expression the properties of the 4-dimensional GB invariant have been 
used (see gT], 02]). 

By varying with respect to the scalar field gives us the equation of motion 



1 -A [V=g (mitW - 2R^F l (<p)d u <p + 7^0)] + ^+ 



V~9 

rlW d P 

— - {Rd^4> - 2R^<f>d»<f>) - -}g = 



(2.7) 



Let us consider the spatially-flat Friedmann-Robertson- Walker (FRW) metric, 

ds 2 = -dt 2 + a{t) 2 (dr 2 + rW) (2. 



where a(t) is the scale factor. Replacing this metric in Eqs. (2.2)-(2.8) we obtain the 
set of equations describing the dynamical evolution of the FRW background and the 
scalar field in the present model: 

H 2 = ^p eff = ^ Q 7 2 + V{4>) + 9i? 2 F a (0)0 2 + 24tf 3 ^) (2.9) 
-2H- 3H 2 = K 2 p eff = k 2 ^ 7 2 - V{(f>) - (3H 2 + 2H^j F 1 

- 2H ( 2F,(0)# + ^) - SH 2 ^ ~ ZH 2 ^ 2 ~ 16 ™^ " 1QH ^ 

(2.10) 



dV , nr T 2 for , dF ll» 



70 + 3 7 #0 + — + 3H Z ( 2F x {<j>)<j> + — ^ + 18^^(0)0+ 

«0 V d(p J ^ n ^ 

12MF X (0)0 + 24 [HH 2 + H^^ = 



In general the couplings -Fi(0) and i*2(0) could be arbitrary functions of the scalar 



field, which gives more general character to the model (2.1) where the couplings 
should be constrained by known observational limits. General couplings also allow 
to increase the number of phenomenologically viable solutions to the DE problem. 
Based on effective limits of fundamental theories like super gravity or string theory, 
the kinetic and GB couplings appear as exponentials of the scalar field (in leading 
a' correction in the case of string theory), but if take into consideration higher order 
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corrections in ol expansion in the effective string theory, the couplings should change. 
Here we will consider the string inspired model with exponential couplings of the form 



F x {(j)) = £ e aK<t,/V2 , F 2 {<f>) = r]e aK<t,/V2 , V{(f>) = y^-^l^ 



(2.12) 



where we also consider an exponential potential that allows to study power-law solu- 
tions. 

In this section we will consider only the standard scalar field corresponding to 7 = 1. 
Lets propose the solution 

H=~, (J) = (J) n \n- (2.13) 



t h 
for quintessential power-law (QPL) expansion, and 

U-t 



H 



t*-r 



) \n 



(2.14) 



for phantom power-law expansion (PPL). By replacing (2.13) and (2.14) in the Fried- 



mann equation (2.9), one obtains the restrictions 

2 <\/2 



an 



3p 2 



1 



'-Ml + v tl + 



9£p 2 2 , A8r]p 3 



'1 



± 



'1 



and from the equation of motion (2.11) we find 



7 (±3p - 1)4> 2 Q - 2V tl + 



6Cp 2 (±3p - 2) 2 48r] P 3 (±p - 1 



'1 



± 



' 1 







(2.15) 
(2.16) 

(2.17) 



where the lower minus sign follows for the phantom solution. Let's consider the two 
cases separately. 

Quintessential power-law expansion 



From Eq. (2.16) and ( J2.1T ) for the quintessence expansion one finds Vq as 

(6£p 2 (3p - 1) + 1 (5p - l)tl)K 2 (j)l + Qp 2 (p - l)t\ 



Vn 



(2.18) 



2(p+ l)K 2 tf 

For 7 = 1, the conditions £ > and p > 1 are enough to give positive potential, i.e. 
Vq > 0. Thus in the regime of accelerated expansion, Vq > 0. One might consider 
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negative kinetic coupling £ < 0, which leads to the following condition by demanding 
positivity of the potential 

6p 2 - l)t\ 



.2i2 



< 



;i - h P )t\ - 6&2(3 P - 1; 

Assuming £ = — 6t 2 (6 > 0), this condition reduces to 

6p 2 (p — 1) 



(2.19) 



.2/2 



< 



(2.20) 



1 — 5p + 66p 2 (3p — 1) 

Taking for instance b = 1, it is clear that in the regime of accelerated expansion 
(p > 1) this condition can be satisfied. 

On the other hand, considering the phantom model (7 = —1), then in the regime of 
accelerated expansion (p > 1) the condition 



K 2 (g < 



6p 2 (p - l)t\ 



0<£< 



(bp 



6p 2 (3p 



or, £<0 



(2.21) 



(5p- l)t\ -6£p 2 (3p- 1 
allows quintessence PL. 

Let's consider only the effect of the kinetic coupling (the effect of the GB coupling 



has been considered in |42j). Setting 77 = in Eqs. (2.16) and (2.17) and solving with 
respect to V and £ one finds 

7(6p- l)K 2 (p 2 + 6p 2 (3p-2) (2p- 7 /t 2 2 )t 2 



V ° 2(3 P +1)kH 2 
A reasonable assumption for 0o could be k 

18p 3 



2p(3p+ l)/t 2 2 



(2.22) 



2^2 



V 



12p 2 + 7 (6p-l) M2 



1 which leads to 
2p-7 



(2.23) 



2(3p+l)t 2 s 2p(3p+l) 1 

In the region of interest that is p > 1 (for both QPL and PPL), we can see that 
independently of 7 the potential is positive. Therefore we have QPL solution for the 
case of pure kinetic coupling (77 = 0). 

A closely related to string theory is the case where V(<p) = 0. To cancel the potential, 



it follows from (2.18) that 



6p 2 (p- l)t 2 



Q& 2 (3p- 1) + 7 (5p- l)tj 



M 2 



(2.24) 



7 



2 >o 



Considering 7 = 1, we note that if £ > then the sufficient condition to have 
is that 1/3 < p < 1. In this case in absence of potential we can not have accelerated 
expansion. More interesting is the case when £ < 0. Let's represent £ = — ht\, where 
6 > 0. Then the condition to cancel the potential becomes 

,2 6P 2 (P - 1) 



—M 2 
I p 



(2.25) 



Qbp 2 (3p — 1) — 5p 

If we take for instance 6=1, then there are two regions of the p-parameter line that 
satisfy the requirement 0q > ; 0.182 < p < 0.633 which is applicable to early-time 
cosmology, and p > 1 which gives accelerated expansion. Then in absence of potential 
it is possible to describe late time QPL expansion in the frame of the present model. 



Considering 7 = —1, if £ < 0, then from (2.24) follows that for p > 1, </>jj > 0, and if 



£ > 0, then the consistency of (2.24) demands that p > 1 and £ < 



(5p-l)tf 
6p 2 (3p-l) - 



In the 



particular case of the model (2.1) without GB coupling, the condition to cancel the 
potential leads to 

'"-'■'-■■^ , _ 7(1 ~ 3p)g (22()) 



6p 2 (2 - 3p) 2 
7(6p-l) p ' 



6p 2 (3p-2) 

in the case 7 = 1, the first equality has sense only for 1/6 < p < 2/3 (see eq. (3.8) 
in [22]), which is appropriate for early-time cosmology. In the case 7 = — 1, the first 
equality is consistent for p > 2/3, which allows accelerated expansion for p > 1. 



Phantom power-law expansion 

For the case of PPL expansion one finds from ( 2.16[ ) and (2.17) (lower sign) for Vo 

(6£p 2 (3p + 1) + 7 (5p + l)tf)K 2 (t)l + 6p 2 (p + l)t\ 



(2.27) 



2(p - 1)kH\ 

For 7 = 1, the potential is always positive for £ > and p > 1. For negative kinetic 
coupling, assuming for instance £ = — 6t 2 (6 > 0) one finds the condition for positive 
potential 

6p2{p + 1) (2.28) 



k 2 4>1 < 



Qbp 2 (3p + 1) — hp — 1 
taking 6=1 for instance, this restriction is consistent for any p > 1. 
Assuming 7 = — 1, then for £ > the potential is positive provided that p > 1 and 
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(setting f = bt\, b > 0) 



QP 2 {P + 1; 



5p — 6bp 2 (3p + 1) + 1 
For £ < and p > 1 (setting £ = — 6t 2 ) one finds 

6p 2 (p + 1) 



(2.29) 



K 



2i2 



< 



Qbp 2 (3p + 1) + 5p+ 1 
Considering only the effect of the kinetic coupling (i.e. rj = 0) one finds 
6p 2 (3p + 2) + 7 (6p+ 1)k 2 0o > (2p + 7 /t 2 



(2.30) 



In 



2p(3p - 1)k 2 (J) 2 



(2.31) 



2(3p - 1)kHI 

For 7 = 1, this potential is positive for p > 1/3 and in this case the kinetic coupling 
becomes negative (£ < 0). Hence, in the case 77 = the PPL expansion takes place 
for negative kinetic coupling. For 7 = —1, the potential is positive for p > 1/3 and 



K 



2a2 



> 



6p 2 (3p+2) 



" 6p+l 

It is interesting also to study the conditions to cancel the potential in the case of 



PPL. From (2.27) follows 



6p 2 (p + l)t 2 



-Mi 



(2.32) 



^ u 7(5p + l)t 2 + 6£p 2 (3p + l)^ p 

For 7 = 1, it follows that for £ > there is not way to cancel the potential for p > 0. 
But if we consider the negative coupling £ < 0, then (setting £ = — bt 2 , b > 0) the 
condition to cancel the potential leads to 

, 2 Qp 2 (p + 1) 



—M 2 
I p 



(2.33) 



^ u 6bp 2 (3p+ 1) - 5p 

taking for instance b — 1, it follows that ( |2.33 ) is consistent for p > 0.483. Therefore 
it is possible to have PPL in absence of potential. Taking 7 = —1 in (2.32), then for 
£ < there is always possible to have PPL for any p, and for £ > 0, the restriction 
(2.32) is consistent provided that £ < 6 ^f^|^ . 

If we limit the model and consider only the effect of the kinetic coupling (rj = 0), 
then from (2.31) for 7 = 1 follows that there is not way to have Vq = 0. Hence, in 
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order to have PPL expansion in this case, it is necessary to have a potential. In the 
case of 7 = — 1, the condition to cancel the potential leads to 

6p 2 (3p + 2)^ 2 & (3p- 



% = M 2 



6p + 1 



p ■ 



l)t\ 



6p 2 (3p + 2) 



(2.34) 



which is consistent for any p > 0. Another important solution of the model (2.1 ) with 



the potential and couplings as given by (2.12) is the de Sitter solution. In fact, if we 
consider the scalar field and the Hubble parameter as constants, i.e. <fi = const. = c 
and H = const. = Hi, then by replacing in (2.9) and (2.11 ) we find 



Hi 



-2c/0 o 



Note that in this solution the kinetic coupling is irrelevant since 
is possible for negative n). 



(2.35) 
(the solution 



In the important the case of V((j)) = 0, as follows from (2.25) and (2.33) there is an 
asymptotic de Sitter solution at p — > oo where 

,o 1 



— M 2 
3b p 



(2.36) 



3 Possible mechanisms to avoid the BR singularity 

The PPL solution suffers the well known problem of the future BR singularity at 
t = t s (133 HH H2J US]. We may use the fact that in the frame of the present model, 
the PPL can be obtained in the absence of potential, and that in the case of dominance 
of potential over the other interaction terms the model presents asymptotic quintom 
behavior. This fact could provide a mechanism to evade the future BR singularity 
as follows. Let's focus on the PPL in the case when we can neglect the effect of the 
potential. To this end we propose the following model 

Fi(0) = £e 2< ^°, F 2 (<t>) = r/e 2 ^ , V(<f>) = V e-^ 1+sw ^\ (6 > 0) (3.1) 

It is clear that the power-law is not a solution to this model, unless 5 = 0. Never- 
theless, we can make some qualitative analysis based on asymptotic behavior of the 
model. In this case, when the curvature is small we assume that the solution behaves 
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as (2.14), which in absence of potential leads to the condition (2.33), that follows for 
PPL expansion. Small curvature means that (t s — t) is large, and the potential that 
behaves like V oc l/(t s — t) 2( - 1+s \ can be neglected compared to the kinetic and GB 
couplings that behave as 1/ (t s — t) 2 . But as the universe evolves, the difference (t s — t) 
becomes smaller and the curvature increases as t — > t s . At this stage the potential 
becomes dominant over the GB and kinetic couplings and the model becomes dom- 
inated by pure quintessential scalar field. In this case the only possible power-law 



solution for the potential (3.1) is of the form 



H 



P 
f 



'° -In* 



1 + 5 U 



which leads to the known conditions 

$ = 2p(l + 5)W, 



V 



P(3p ~ 1) M 2 



t 2 
1 1 



(3.2) 



(3.3) 



giving an EoS parameter w > — 1, avoiding in this way the future BR singularity. 
As the universe continue evolving after the dominance of the potential, the curva- 
ture turns again to small values and the interacting terms start dominating again. 
Nevertheless, when the potential and the couplings are present there is an important 



solution, namely the de Sitter solution. In fact, if we assume for the model (J3 . 1 p 

<p = const. = c, H — const. = Hi 



then, by replacing in (2.9) and (2.11) this solution gives 

1 



3 -2c/?Jo 



In 



8r]K 4 V 



(3.4) 



(3.5) 



8nK 2 ' 25 

This is an alternative to the approach presented in jl2] where the GB coupling and 
phantom scalar field were considered. If we consider the phantom version of the model 
(i.e. with 7 = — 1), then there are more alternatives to avoid the BR singularity. 



Phantom scalar field 7 = — 1 



Bellow we consider the phantom kinetic term in the model (2.1) (i.e. assuming 
7 = —1), which gives more possible ways of evading the BR singularity, namely 
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1. In absence of kinetic coupling = 0). 
In this case, the model reproduces the same results presented in |42j. 



2. In absence of GB coupling (rj = 0). 
One may consider the situation when initially at low curvature (large time) the po- 
tential term dominates during PPL expansion, and then when the solution is closer 
to the BR singularity the kinetic coupling becomes dominant, allowing the possibility 
of QPL expansion with EoS parameter w > — 1. To this end we propose the model 

V = Voe- 2 ^ , Fx = £e 2(1+5) ^°, (6 < 0) (3.6) 



Assuming that the solution behaves as (2.14), then neglecting the kinetic coupling, 



from (2.9) and (2.11) (for 7 = —1) follow the known conditions for existence of PPL 



solution (2.14) 



% = 2pM 2 , V = P{3P + 1] Ml 



as the curvature increases when t — > t s , the kinetic coupling becomes dominating 
and the potential could be neglected. In this case we may assume the same solution 



(3.2) for the quintessential expansion, that being replaced in (2.9) and (2.11) (setting 

3p — 



V = and F 2 = 0) leads to the conditions 
. 6p 2 (3p-2) 



(i + 5)X, f 



-ti 



(3i 



6p-l "~ ' ~' "' pi " 6p 2 (3p-2) 
which is consistent for p > 2/3. So there is possible to change the effective EoS 
from w < — 1 to M) > - 1 avoiding the BR singularity. Hence, when the term with 
kinetic coupling becomes dominant, the BR singularity might be prevented. So the 
kinetic coupling may play a role similar to the GB coupling [42] in working against 
the singularity. 



For the power-law solution of the form (3.2), after the dominance of the kinetic 
coupling the curvature becomes small again, and the potential term recovers his 
dominance, but there exists a solution when both, the potential and the kinetic 
coupling are present which corresponds to a de Sitter phase. If we assume the solution 



(3.4), then replacing in (2.9) and (2.11) with (3.6), one finds 



Hi 



K 2 V 



-2c/4> 



(3.9) 
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which is valid for arbitrary c, since = 0, making this solution independent of the 
kinetic coupling. Therefore, there is a possibility that the universe enters in a de 
Sitter phase after domination of the kinetic coupling. 



3. In absence of potential (V = 0). 
Here we consider the following model 



Fi = £e 



20/0 o 



rje 



2(l+5)0/0o § < o 



(3.10) 



Assuming the phantom solution (2.14), we see that at low curvature when t s — t is 



large, the GB coupling behaves as l/(t s — t) 2+2S and can be neglected with respect 



to the term with kinetic coupling. In this case, by solving Eqs. (2.9) and (2.11) one 
finds the conditions 



6p 2 (3p + 2) 



Ml 



3p + l _ t2 



(3-11) 



Qp + 1 "~ pi " Qp 2 (3p + 2) 1 

which corresponds to phantom phase with effective EoS w < — 1. When the curvature 
turns to large values at t — > t a , the term with kinetic coupling could be neglected 
giving rise to the dominance of the GB term. Neglecting the kinetic coupling in 



(3.10) and assuming a solution of the form given by eq. (3.2), one finds from (2.9) 



and (2.11) 



6p 2 (p — 1) 



1 + 5) 2 M 2 



v 



V 



8p(5p-l) pl 



(3-12) 



5p-l 

which is consistent for p > 1, leading to quintessential expansion with w > — 1. 
Note that the role of the kinetic and GB couplings could be changed, i.e. one starts 
with domination of the GB term and ends with domination of the kinetic coupling, 
by proposing 



F 1 = £ e 2(1+5)< ^°, F 2 =r/e 20/00 , 5 < (3.13) 
In this case, at low curvature the kinetic coupling can be neglected, and replacing the 



3 20/0o 



PPL solution (|2.14|) in (|2.9|) and (|2.11|) one finds the conditions 

3p + l 



2 - 6p2(p + 1) Ml V 



M 2 t 2 



5p+ 1 



8p(5p + 



(3.14) 
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which is consistent for any p > and negative r\. At large curvature when t — > t s , the 



kinetic term becomes dominant and we can propose the solution (3.2), which being 



replaced in (2.9) and (2.11) gives the restrictions 



6p 2 (3p- 2) 



3p- 1 



(3.15) 



6p-l K ' • '"P 7 * 6p 2 (3p-2) 

which is consistent for accelerated expansion with p > 1 and effective EoS io > — 1. 
However in the model without potential, there is not de Sitter solution corresponding 
to constant scalar field. 



4- With all the terms. 
We may consider the phantom scalar model with the potential and couplings given 
by 

F 1 ( ( j ) ) = £e 2(1+5) ^°, F 2 {<j>) =r/e 2(1+5) ^°, V{<j>) = V^ 1 ^, (5 < 0) (3.16) 
If we neglect the couplings Fi and F 2 at low curvature, then the only possible power- 



law solution is the phantom one, given by (2.14). When the curvature becomes large 



at t — > t s , the interacting terms become relevant and (neglecting the potential) there 



is a power-law solution of the form (3.2) which leads to 

m2 6p 2 (p-l)tj 



(1 + 5YM 2 p 



(3.17) 



ru (5p -l)t\- 6p 2 £{3p - 1) 

which is positive whenever t\ > 6£p 2 (3p — l)/(5p — 1). For negative £ this condition 
is satisfied for any p > 1, leading to w > — 1 and avoiding the BR singularity. Note 



that the solution (3.17) exists in the asymptotic de Sitter limit at p — » oo («?—»■ — 1), 
given by 

-M 2 V (3.18) 



3£ 



valid for £ < (in this limit 77 — > 0). After the dominance of the coupling terms the 



curvature begins to decrease again, but there exists a de Sitter solution (3.4) when 



the three terms in (3.16) are present 

1 



2(l+5)c/4> 



5 7/ fir 



00 

25 



In 



8t7k 4 1/ 



(3.19) 
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So it is possible to implement the asymptotic mechanism to avoid the BR singularity 
as proposed in [42J, in different variants of the model depending on the correlation 
between the kinetic coupling, the GB coupling and the potential. As has been shown, 
this mechanism can be implemented in the standard and phantom version of the 
scalar field. 

It is worth mentioning that the account of quantum effects near the singularity were 
also considered to moderate the BR singularity [47]. Another interesting alternative 
to avoid the future BR singularity is provided by the solutions known as "Little 
Rip" (LR) |3SJ HHJ [50j 36J, which are free of future singularity. The LR solutions 
produce late-time cosmological effects similar to that of the BR solutions, as the 
rapid expansion in the near future with an EoS w < — 1, but the scale factor and 
density remain finite in finite time. As in the case of BR, the LR solutions also lead 
to the dissolution of all bound structures in the universe in the future. 



4 Stability of the power-law solution for the string 
motivated model 

We use the dynamical system approach in order to analyze the stability of the above 
power-law solutions, in the specific case of V — 0, that apart from simplifying the 
dynamical system is also closely related to string theory. Let's introduce the following 
dimensionless variables: 

K(j) 2 '2 

x = -=- : , k = 3k F x (p , 

V2H . (4.1) 

2 -dF 2 H 

9 = 8KH(f) W € = 1p 

In fact, these variables are related with the density parameters for the different sectors 
of the model 



where 



1 dF 

-(0 2 + n Pk = §K 2 H 2 F^\ Pgb = 2Ak 2 H 3 ^ (4.3) 
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The Eq. (2.9) imply the following restriction on the density parameters 



GB 



(4.4) 



and the effective equation of state (EoS) can be written as 

W e ff — W^VL^ + W k ^lk + Wgb^GB = "I 



(4.5) 



Introducing the e-folding variable iV = log a, in terms of the variables (3.1), the Eqs. 



(2.9)-(2.12) can be transformed into the following first-order autonomous system 



7X 2 + 3k + 3g - 3 = 



2-fxx' + 2 7 (3 + e)x 2 + k' + 2(3 + 2e)k + 3(1 + e)g = 
2e + 3 + 7a; 2 - \k' - -(3 + 2e)k - g' - (2 + e)g = 



k' = (ax + 2e)k + 2—k 



x 



x 



9 



(ax + 2e)g H g 

x 



(4.6) 

(4.7) 
(4.8) 

(4.9) 

(4.10) 



where '"" denotes derivative with respect to N and 7 = ±1 is the sign of the free 
kinetic term. Note that the last three Eqs. come from the explicit form of the 



potential and the couplings given in (2.12). From Eq. (4.8) follows the expression for 
the slow-roll parameter e 



9 + 3 7 x 2 - 2k' -3k- 3g' - 6g 
2k + 3g-6 



(4.11) 



It is easy to check that the power-law solutions (2.13) and (2.14) are critical points of 



the system, i.e., if we write the dynamical variables (|4.1) for H and 4> given by (2.13) 

16r]K 2 p 



or (2.14) as 



K(p , 3£k 2 4> 2 



k 



t 2 
1 1 



9o 



t 2 
'1 



(4.12) 



where the "-" sign is for the PPL, then these variables satisfy the equations: x' 
^0 = 9o = 0- S° we will consider small perturbations 



x = Xq + 5x, k = k + 8k, g = go + Sg 



(4.13) 
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and check the stability around the critical point (xo, k , g Q ). Solving the system (4.6)- 
(4.11) with respect to x', k',g' one can write 



x 



fi{x,k,g), k' = f 2 {x,k,g), g' = f 3 (x,k,g) 



(4.14) 



for small perturbations, x', k',g' suffer the change 

/ sx' \ /^A f 
5k' 



\5g> J 



dfl 
dx 

\ dx 



9h 

dk 

dfl 
dk 

dh 

dk 



9g 
dh 
9g 

dh 
dg 



J 



( 8x \ 

5k 
\Sg J 



(4.15) 



where the matrix is valuated at the fixed point (xq, kg, go) given by (4.12). The sta- 



bility under small perturbations demand that the eigenvalues of the above matrix be 
negative or complex with negative real component. We will analyze the stability for 
two different cases. In the first case we consider the model with only kinetic coupling 
(g = 0), and in the second case we consider both couplings (the case with only GB 
coupling was considered in [42J). 



The model with non- minimally coupled kinetic term. 



By setting g = in (4.6)-(4.11 ), it reduces to a two dimensional system for x and fc, 



and for small perturbations we may write 





(4.16) 



where the two dimensional matrix M is evaluated at the critical point (xq, ko), which 
gives the components 

3(5^0 T 9)^o ~~ (6ax - 9)k^ + (±72xq - Qax =F Sax^ + 9)k 



M 



12 



6fcg ± 6x 2 - 2(±x 2 - 3)k 
(±24xq — 3ax =F 2ax^ — 6(ax — 3)k + 9)x 



M 



21 



6fcg ± 6x 2 - 2(±x 2 - 3)k 
3((axo — S)ko — 3axo + 12)/c a;o 



3/lq ztz 3Xq (^i^Q 



3)fc 



M> 



22 



(24/c + 3axQ — 2akoXo — l&)x\ 
3k 2 ±3x 2 - (±x 2 -3)k 



(4.17) 
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where the lower sign is assigned to the phantom model (7 = —1, see Eqs. (2.9) and 



(2.11)). Replacing x an d k from (4.12) and using the restrictions (2.26) for QPL 



we find the following eigenvalues for the matrix M 



Ai 



3p 



V 



A 2 



3p 



V 



- < p < - 
6^3 



(4.18) 



where p should obey the above restriction for consistency, according to (2.26). In 



the interval 1/6 < p < 2/3, A2 > and therefore the power-law solution of the form 
H = p/t is unstable for the model in absence of potential and GB term. Note that the 
PPL solution with only kinetic coupling is not possible as was demonstrated above 
(see eq. fl2~3lj) for V = 0). 



On the other hand, if we consider the phantom model that obeys the Eqs. (2.9) and 



(2.11) (7 = —1), then taking into account the lower sign in the components (4.17) 



one finds the following eigenvalues for the QPL 



Ai 



l-3p 

? 

p 



A, 



-3p 2 
p 3 



(4.19) 



where the last inequality follows from the consistency of the solution of Eqs. (2.9) 
and (2.11) with V = and 77 = 0, for QPL (see eq. (3.8)). Therefore, for the 



phantom model with kinetic coupling, the power-law solution ( |2.13 ) is a stable fixed 
point provided p > 2/3. 



If we consider the phantom model with PPL (2.14), then using the lower sign in (4.17) 



and in (4.12) one finds the following eigenvalues 

3p + 2 



Ai 



A, 



V 



3p+l 
p 



(4.20) 



which is valid for any p > 0, as follows from Eqs. (2.9) and (2.11) with V = and 
77 = 0, for PPL. Then the power-law solution (2.14) is a stable fixed point for the 
phantom model with V = and 77 = 0. 



The model with kinetic and GB couplings 

Here we consider the stability of power-law solution for the three dimensional au- 



tonomous system (4.6)-(4.11) in the following cases: 
Quintessence and phantom power-law for 7 = 1. 
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Evaluating the matrix elements of eq. (4.15) for the fixed point (4.12) (upper sign) 
and for the QPL (upper sign), we find the following eigenvalues 

1 - 3p . 2-3p 



A 1 = 0, A 2 



A s 



(4.21) 



p p 

note that A 2 and A3 are negative for p > 2/3. This restriction is compatible with the 



condition of consistency for the QPL solution (2.13), as can be seen in (2.25). 
For the PPL we have found the eigenvalues 

2 + 3p x 1 + 3p 



Ai = 0, A 2 



A, 



(4.22) 



p p 

where A 2 and A3 are negative for any p > 0, which is compatible with the restrictions 
Quintessence and phantom power-law for 7 : 



■1. 



Evaluating the matrix elements in (4.15) for the fixed point (4.12) (lower sign) and 
for the QPL (upper sign), we find the following eigenvalues 

1 - 3p . 2 - 3p 



Ai = 0, A 2 



A, 



V 



P 



(4.23) 



which is compatible with the solutions in absence of potential for p > 1 (see 3.17). 
For PPL we have the eigenvalues 

2 + 3p . l + 3p 



A 1 = 0, A 2 



A s 



(4.24) 



p p 

valid for £ < and p > 0. Note that the eigenvalues are independent of 7. So in the 
presence of both couplings the stability of the power-law solution does not depend 
on the standard or phantom character of the model. But a problem appears in both 
cases due to the presence of zero eigenvalues. In this case the linear expansion fails 
to provide information on the stability of the fixed point. We need to consider higher 
order corrections to study the stability of perturbations along the zero eigenvalue 
direction. 



The centre manifold analysis. 

To analyze the stability in the presence of zero eigenvalues we use the approach of 
the central manifold |51j . [52] , [53J, which reduces the dimensionality of the system 
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near the critical point, and limits the stability analysis to the reduced system. Thus 
the stability properties of the system become determined by the (in)stability of the 



reduced system. To this end we need to translate the fixed point (4.12) to the origin, 
by introducing the variables (we keep the same symbols) 



x ->■ x - x , k -» k- k , g ->■ g - g 



(4.25) 



The simplest case takes place for one zero eigenvalue, which leads to one-dimensional 
reduced system. Composing the matrix M$ with the eigenvectors of the Hessian in 
the new defined fixed point (0,0,0), we introduce a new set of coordinates 

I u\ I x \ 



V 



M 



k 

\ 9 ) 



(4.26) 



In these coordinates the dynamical equations can be written in the form 

' \ f u\ ( h(u,v,w) \ 



1 


v! 


\ 






v' 






\ 


w' 


/ 













l-3p 
P 

2-3p 



V 

\w ) 



+ 



f 2 (u,v,w) 
V h(u,v,w) J 



(4.27) 



where the last column represents the non-linear terms. Note that in our case the 
variable u is actually the same variable x. Then the system can be written as 



u' = fi(u,y), 
V = Ly + f(u, y) 



where 





f V \ 




f l-3p 


H 






V 


K W j 




^ 





2-3p 



f{u,y) 



h{u,y) 
h(u,y) 



(4.28) 
(4.29) 



We now turn to the definition of centre manifold: 
The space 



W c {0) = {{u,y) e R 1 x R 2 \y = h(u), \u\ < S,h{0) = 0,Dh{0) = 0} 



(4.30) 



where Dh is the matrix of first derivatives of the vector valued function y = h(x), is 



called the centre manifold for the system (4.28). Since y = h(x), the dynamics of the 
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system becomes reduced to the centre manifold in the neighborhood of x, and the 
stability properties of the full dynamical system depend on the analysis in the centre 



manifold. Using y = h(x), the system (4.28) leads to 



Dh{u)h{u, h{u)) = Lh{u) + f{u, h{u)) 



(4.31) 



This differential equation can be used to find h(u), and then by replacing h(u) into 
the first equation (4.28) (i.e. v! = fi(u, h{u))) we can analyze the stability of the 



reduced system. Near the critical point we can Taylor expand h(u) in powers of u 



and calculate the coefficients of the first non-trivial terms from (4.31). In our case we 
assume h of the form 



a 2 u 2 + a 3 u 3 + a^u 4 + a 5 u 5 + 0(u 6 ; 
h{u) = (4.32) 
V b 2 u 2 + b 3 u 3 + b 4 u 4 + b 5 u 5 + 0(u 6 ) J 

Using the restrictions on £ and r\ for the QPL solution in absence of potential, the 



critical point (4.12) takes the form 

3«V - 3«V + 20p - 4 



x 



2 



a 2 (3p — l)p 2 



9o 



2(9a V - Qa 2 p 2 + 24p - 4) 
3a 2 (3p — l)p 2 



(4.33) 

where a is given by (2.15). To avoid large analytical expressions we will limit the 



analysis to specific values of p. The first interesting value of p corresponds to the 
de Sitter limit, which gives the critical point (x = 0,k = —l,go = 2) and the 
corresponding eigenvalues from (4.21) are (Ai = 0,A2 = — 3, A 3 = —3). Applying the 



above central manifold analysis, and after large but straightforward calculations we 
find the following equation for the reduced system 

9 3 
64 



u 



16a 2 



u 7 + 0(u* 



(4.34) 



by replacing h(a) 



16a 2 



Jj and integrating this equation one obtains 



u = u Q (l - 6u 6 Q h(a)N) 



-1/6 



(4.35) 



where uq is the initial perturbation along the zero eigenvalue direction. In order for 



this initial perturbation to decay it follows from (4.35) that independently of the sign 
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of uq, the coefficient h(a) must be negative. Therefore if h(a) < the critical point 
will be stable. Turning to eq. (4.34) we see that if a 2 > 12 then h < and the critical 



point is a stable attractor. 

Another reasonable value is p = 200/9, which gives the EoS parameter w = —0.97. In 
this case the critical point is (x = k = -|g - 19 ^ 7 0a2 , 9o = f§? + qssooo^ ) and 
the eigenvalues are (A x = 0, A 2 = —591/200, A 3 = —291/100). The centre manifold 
analysis gives 

u = fi(a)u 7 + 0(u s ) (4.36) 

but the analytical expression for p(a) is too large to be displayed here. We limit 
ourselves to numerical evaluation for some values of a: (a = 2, /i = 0.809), (a = 3, 
fi = 0.03526), (a = 12, fi = -0.03592) and (a = 14, fi = -0.03549). Note that for 
the last two values of a the critical point is stable. 



One can also perform the same analysis for the PPL (lower sing in (4.12) with eigen 



values given by (4.22). For the specific value of p — 80/3 giving the EoS parameter 



w = —1.025, the centre manifold analysis gives 

u = u{a)u 7 + 0(u 8 ) (4.37) 

evaluating u(a) for some values of a, gives: (a = 1, v — 0.055), (a = 3, v — —0.008), 
(a = 12, v — —0.049), (a = 14, v = —0.051). The last three values of a lead to 
stability. 



5 Discussion 



We studied late time power-law cosmological solutions based on string spired scalar- 
tensor model including a coupling to the Gauss-Bonnet invariant and kinetic couplings 
to curvature. The model allows quintessential and phantom power-law expansion in 
a variety of scenarios that involve different asymptotic limits. In the case with poten- 



tial the conditions (2.18) (2.19) and (2.21) allow QPL expansion, and the conditions 



(2.27), (2.28), (2.29) and (2.30) allow PPL expansion. In absence of the GB coupling, 



the restrictions (2.22) allow QPL and the restrictions (2.31) allow PPL. When we 
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neglect the potential (which is closely related to string theory), the condition (2.24) 
allows QPL expansion and ( |2.32 ) allows PPL. If in addition to the potential we ne- 
glect the GB term, then the model continues to have power-law solutions, where in 
this case for 7 = 1, the condition (2.26) allows power-law but in the range of decel- 



erated expansion. And for the case of 7 = —1, it follows from (2.26) that the model 
has QPL solution. Concerning PPL solutions in this limit (i.e. V = 0, 77 = 0), then 
from (2.31) for 7 = 1 it follows that there is not PPL without potential, and in the 



case of 7 = —1, the condition (2.34) allows PPL for any p > 0. 

The model also exhibit de Sitter solution in various scenarios: considering = const. 



and H = const., then from (2.9) and (2.11) follows the solution (2.35) and in absence 



of potential the asymptotic de Sitter solution (2.36) takes place. We also investigated 
the possible mechanism to avoid the Big Rip singularity in the case of PPL expansion. 
To this end we made a qualitative analysis, by proposing different asymptotic scenar- 



ios where one or two interacting terms (including the potential) in (2.1) are relevant 
at low curvature, while the remaining terms become relevant at large curvature (when 
t — > t s ), providing a quintessential solution and evading in this way the singularity 
(see [42J). After that, the universe might evolve asymptotically towards a de Sitter 
solution as was shown in different scenarios. The mechanism was implemented in the 
standard and phantom version of the scalar field (i.e. 7 = ±1). Of special interest is 
the first case (3.1) (7 = 1) where at t — > t s the dynamics becomes dominated by pure 
quintessential scalar field, which is the only possibility. 

We have performed the stability analysis for the string inspired model (i.e. in absence 
of potential) and have found that the power-law solution (2.13) (or ( 2.14[ )) is a critical 
point of the model and is stable fixed point in different scenarios: we first considered 
the case where the GB coupling is neglected (i.e. V = 0, r\ — 0). In this case, for 
7 = 1, the critical point is unstable in the allowable range of p (1/6 < p < 2/3, 



see (2.26)). For 7 = —1, the critical point (2.13) is stable attractor, provided that 
p > 2/3. The PPL solution (2.14 is a stable attractor for the phantom model (7 = — 1) 
for any p > 0. The scalar field with GB correction (i.e. with V = 0, £ = 0) was 



considered in (42]. Next we analyzed the stability of the model with GB and kinetic 
coupling terms. In the model with 7 = 1, for QPL we have found the eigenvalues 
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(0, (l-3p)/p, (2-3p)/p) and for PPL (0, -(2 + 3p)/p, -(1 + 3p)/p). For 7 = -1, we 
have found the same eigenvalues. Due to the presence of zero eigenvalues, the linear 
expansion fails to provide information on the stability of the fixed point. We used 
the centre manifold analysis to determine the stability properties of the critical point, 
and found that in the important limit of de Sitter solution the critical point is a sta- 



ble attractor, under certain restriction coming from the expression (4.34). Numerical 
evaluation was also done for concrete values of p and it was found stability for some 
cases. 

Resuming, the presence of the kinetic coupling besides the GB coupling, extends the 
number of possible scenarios to realize cosmological solutions with BR singularity, 
compared to the model with only GB correction. Additionally, the different asymp- 
totic cases of the present model not only extend such possibilities, but also provide a 
number of alternatives to avoid the BR singularity leading to an universe that might 
evolve towards a de Sitter phase. Of special interest is the first case given by the 



model (3.1 ), in which the BR solution is obtained without appealing to phantom field 
(i.e. 7 = 1), and near the singularity the dynamics becomes dominated by purely 
quintessential scalar field with EoS w > — 1. Note that in this scenario the terms with 
couplings produce the PPL and the scalar potential acts against the BR singularity, 



while in the other cases (7 = —1, see Eqs. (3.6), (3.10), (3.16)) the GB and kinetic 
couplings might prevent the BR singularity. The above results show that the string 
effects could play significant role in late time cosmology. 
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